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I ; ABSTRACT 

^ ; 

Oh' Context. Our previous models of a giant planet migrating through an inner protoplanet/planetesimal disk find that the giant shepherds a portion 
^ of the material it encounters into interior orbits, whilst scattering the rest into external orbits. Scattering tends to dominate, leaving behind 
» I abundant material that can accrete into terrestrial planets. 

Aims. We add to the possible realism of our model by simulating type I migration forces which cause an inward drift, and strong eccentricity 
and inclination damping of protoplanetary bodies. This extra dissipation might be expected to enhance shepherding at the expense of scattering, 
' ' ■ possibly modifying our previous conclusions. 

Methods. We employ an N-body code that is linked to a viscous gas disk algorithm capable of simulating: gas accretion onto the central star; 
gap formation in the vicinity of the giant planet; type II migration of the giant planet; type I migration of protoplanets; and the effect of gas 
drag on planetesimals. We use the code to re-run three scenarios from a previous work where type I migration was not included. 
Results. The additional dissipation introduced by type I migration enhances the inward shepherding of material but does not severely reduce 
\^ scattering. We find that > 50% of the solids disk material still survives the migration in scattered exterior orbits: most of it well placed to 
complete terrestrial planet formation at < 3 AU. The shepherded portion of the disk accretes into hot-Earths, which survive in interior orbits for 
the duration of our simulations. 

Conclusions. Water-rich terrestrial planets can form in the habitable zones of hot- Jupiter systems and hot-Earths and hot-Neptunes may also 
be present. These systems should be targets of future planet search missions. 

o : 

|L* ' Key words, planets and satellites: formation - methods: N-body simulations - astrobiology 
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^ 1. Introduction. Lin et aL 1996; Ward 1997; Nelson et al. 2000). According 

■ - - to recent models of the late stages of solar system formation, 

The core accretion hypothesis of giant planet formation pos- , • ■ . i ,.u * » f *u i 

-'^ & F f which appear to explain the current architecture of the solar 

tulates that giant planets originate in a cool location in a pro- . .rp . . ^ , onnc: a/t u j n- » i iaac 

^ f s= F system (Tsigams et al. 2005; Morbidelh et al. 2005; Gomes et 

toplanetary disk, beyond the snowline where the condensation , oaacn i j c ... » u j i 

^ al. 2005), Jupiter and Saturn appear to have undergone only 

of ices boosts the mass of solid material present (e.g. Pollack et i - •-. j ■ -.• ^- tu ^ <.u u f a 

^ ^ ° limited migration, suggesting that they may have formed quite 

al. 1996; Papaloizou & Nelson 2005). This, and the fact that re- i -.u <.<. a- i t^u t-u • • « 

' ^ ' late with respect to dispersal of the gas. Their influence on inte- 

moteness from the central star enhances the gravitational reach .. » ii »r ^ a- ^ ^tatt u 

" nor terrestrial planet formation at distances < 2 AU was prob- 

of accretion, may allow a ~ 10 core to grow and then ac- , , a ^ uu u » j i » i ■ » 

^ w & ^l,jy modest, although recent models suggest secular interac- 

cumulate a massive gaseous envelope before the loss of the .. , i j • » » i j • a- i a- 

° ^ tions may have played an important role during gas disk dis- 

nebular gas which is observed to occur ~ 1-10 Myr after star , . , onnc\ t » » i 

" ^ persal (Nagasawa et al. 2005). In contrast, many giant exoplan- 

formation (Haisch et al. 2001). The alternative gravitational in- . , f u * i • u • * j n 

^ ^ " ets show evidence for substantial migration, being found well 

stability model also predicts that giant planets will form in the .... , . , ■ , , ... » u i ' r 

^ ^ fe F Within what would have been their parent nebula s snowline 

cool outer disk regions (Boss 2000). Giant planets are also ex- /d »i » i oaa^n »u » * i u • i 

° ^ ^ F (Butler et al. 2006), the most extreme examples being known 

pected to migrate inward from their original formation location u » i -^u i * j » » ^ a i att 

f ^ & as hot-Jupiters with examples found at distances < 0.1 AU 

due to gravitational interactions with the disk, either during the , -. -.f -.n a - i» »u a 

° ' ° down to just a few stellar radii. These planets must have passed 

core stage via type 1 migration, where the body remains embed- u *u • • » * i / i ia\^ n j 

° -'^ ° through their inner systems at an early time (~ 1-10 Myr) and 

ded in the gas (e.g. Ward 1997; Papaloizou & Larwood 2000; , ■ • v ui u u j . a f 

° ^ ° ' F would inevitably have had a strong influence on the formation 

Tanaka et al. 2002; Tanaka & Ward 2004), or via type II migra- ^ i i » »u » u »u • u ^ a^ »u 

' ^ ^ ° of rocky planets that would otherwise be expected to grow there 

tion when the planet has grown massive enough (> 100 Me) ^„ timescales of ~ 10-100 Myr (Chambers 2001 ; Kleine et al. 

to open a gap in the gas and becomes coupled to the inward 

viscous evolution of the disk (e.g. Lin & Papaloizou 1986; 
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2002; Halliday 2004; O'Brien et al. 2006; Niirnno & Kleine 
2007). 

It is expected that terrestrial planets external to a hot- Jupiter 
would have stable orbits (Jones et al. 2005) and should be 
able to form from any available pre -planetary material with a 
period ratio > 3 (Raymond et al. 2005a). However, until re- 
cently it has been a common approach in both astrophysical 
and astrobiological literature to assume that terrestrial planets 
could not have formed in hot- Jupiter systems due to the disrup- 
tive effect of the giant planet's migration which is deemed to 
have cleared the inner system of planet-forming material (e.g. 
Ward & Brownlee 2000; Lineweaver 2001; Armitage 2003; 
Lineweaver et al. 2004). However, recent work by two groups 
who have modeled terrestrial planet formation in the pres- 
ence of type II giant planet migration have shown that this as- 
sumed clearing does not happen (Fogg & Nelson 2005, 2006; 
Raymond et al. 2006; Fogg & Nelson 2007; Mandell et al. 
2007). Their simulations show that the giant planet shepherds 
the solids disk inward, compacting it and exciting the orbits of 
objects caught at first order resonances. Much of this excited 
material eventually experiences a close encounter with the gi- 
ant planet and is expelled into an exterior orbit. The net effect 
of the giant's passage is therefore not the destruction and dis- 
appearance of the inner system disk, but instead a modest dilu- 
tion and strong excitation and radial mixing of solid material. 
Over the entire course of the migration, a new disk of solid 
material builds up in orbits external to the final position of the 
hot- Jupiter where terrestrial planet formation can resume. The 
results of further simulation of accretion in this scattered disk 
suggest that terrestrial planets should be commonplace in hot- 
Jupiter systems: worlds that might possibly be water-rich due 
to strong inward mixing of material from beyond the nebular 
snowline (Raymond et al. 2006; Fogg & Nelson 2007; Mandell 
et al. 2007). 

All these papers agree that the outcome of giant planet mi- 
gration through an inner system soUds disk is to partition most 
of its mass into the exterior scattered remnant described above 
and an interior remnant resulting from shepherded material that 
is not scattered by the giant. This shepherded fraction is com- 
pacted into the restricted volume close to the central star, where 
dynamical times are short, resulting in a rapid burst of accre- 
tion with the potential to form stable hot-Neptune or hot-Earth 
type planets orbiting interior to the final position of the giant 
(see also Zhou et al. 2005). Comparing results however reveals 
that the relative predominance of scattering versus shepherding 
is significantly influenced by the strength of dissipative forces 
(gas drag and dynamical friction) operating in each particular 
model. Whilst interior and exterior remnants are found in all 
cases, a distinct trend is evident favoring shepherding under 
conditions of strong dissipation and scattering when dissipation 
is weaker. For example, the model in Fogg & Nelson (2005) 
assumed a steady-state gas disk with a fixed r~^-^ surface den- 
sity profile which maximized both the influence of gas drag on 
planetesimal sized objects and hence their ability to damp the 
orbits of protoplanets via dynamical friction. The formation of 
shepherded hot-Neptunes and hot-Earths was a common out- 
come of these simulations with as little as ~26% of the orig- 
inal disk material being scattered into the external disk in the 



youngest scenarios. (In the case of constant gas density the av- 
erage strength of dynamical friction still falls with disk matu- 
rity as planetesimals are cleared through accretion.) In contrast, 
in their latest model (Fogg & Nelson 2007) which incorporates 
a time dependent, viscously evolving, gas disk with a partial 
central cavity and annular gap centered on the giant planet's 
orbit, dissipation is much lower and declines with time, espe- 
cially in the inner regions of the system. In this perhaps more 
realistic case, hot-Earths start to form in the shepherded portion 
of the disk but are almost always scattered or accreted by the 
giant at late times when it has migrated down to a < 0.5 AU. 
Irrespective of the maturity of the solids disk at the start of 
migration, > 60% of its mass is always scattered into an ex- 
ternal disk. Dissipation levels within the model of Raymond et 
al. (2006) and Mandell et al. (2007), with its linearly declining, 
fixed surface density profile gas disk, appear to lie intermediate 
between the two models of Fogg & Nelson. Shepherded hot- 
Earths sometimes survive at the end of their simulations and 
sometimes not, although a rigorous comparison between their 
work and that of Fogg & Nelson is not possible due to differ- 
ences in assumed initial conditions, and the fact that in most 
of their simulations the giant planet comes to rest at a radial 
distance of > 0.1 AU. 

One potentially important physical process neglected in 
these models is type I migration which operates on sub-gap 
opening bodies of ~ 0.1 - 100 Me (e.g. Ward 1997; Papaloizou 
& Larwood 2000; Tanaka et al. 2002; Tanaka & Ward 2004). 
Type I migration is generated by the asymmetric torques ex- 
erted on a protoplanet from the wakes it creates in the gas disk 
and is expected to impart an inward radial drift and strong ec- 
centricity and inclination damping that increases in effective- 
ness linearly with protoplanetary mass. The realism of type 
I migration has long been controversial as theoretically pre- 
dicted spiral-in times can be so short as to threaten proto- 
planetary survival and to render an explanation of the solar 
system's architecture particularly difficult. Many planet for- 
mation models therefore ignore this process entirely. It might 
be that type I migration does not operate as rapidly as pre- 
dicted: one proposal being that stochastic torques from density 
fluctuations in a turbulent disk could superimpose a random 
walk in semi-major axis over the smooth inward drift predicted 
by theory (Nelson & Papaloizou 2004; Laughlin et al. 2004; 
Nelson 2005). However, it is becoming apparent that even near- 
nominal rates of type I migration are not necessarily fatal to 
planet formation and survival. Work by McNeil et al. (2005) 
has shown that it is possible to grow a terrestrial planetary sys- 
tem in a simulation including type I migration forces by en- 
hancing the original quantity of solid material present whilst 
including a rapidly dissipating gas disk. The model of Daisaka 
et al. (2006) succeeds in retaining sufficient material in the ter- 
restrial region by invoking a shorter gas depletion timescale 
in the inner system than the observational value (based on in- 
frared observations at a ~ 100 AU), or if the the initial gas to 
dust ratio is smaller than the conventional minimum mass solar 
nebula (MMSN) model. The problem with preventing the more 
massive giant planet cores from being lost to the central star is 
of greater severity, but recent work by Thommes & Murray 
(2006) suggests that this may be possible late in the lifetime of 
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the gas disk when accretion and migration timescales become 
comparable. Proposed system properties that might widen such 
a window for core survival include small planetesimal sizes, 
low midplane gas disk viscosities and enhanced collision cross 
sections due to core atmospheres (Chambers 2006a; Thommes 
et al. 2007). 

The possible realism of strong type I migration forces is of 
obvious and perhaps critical relevance to the above-cited mod- 
els of giant planet migration through an inner system sohds 
disk. Even if the issue of pre-existing type I migration of the 
giant planet's core is neglected and the assumption of the fully 
formed giant planet undergoing type II migration is retained, 
type I migration forces would still be expected to operate on 
the protoplanetary components of the inner system disk, in- 
troducing another source of damping additional to dynamical 
friction. Extra dissipation would be expected to affect the par- 
titioning of inner system material by enhancing the shepherded 
fraction at the expense of the scattered fraction, potentially in- 
vaUdating all the previous claims that sufficient solid material 
survives the passage of the giant to permit the formation of ter- 
restrial planets. Our primary motivation is to investigate this 
issue here by adding type I migration forces to the model of 
Fogg & Nelson (2007) and re -running one early and two late 
scenarios from that paper. We note that additional sources of 
damping may also be important such as damping by coUisional 
debris, but inclusion of these effects is beyond the scope of our 
model. In Section 2 we outline the additions to our model and 
the initial conditions of the simulations; in Section 3 the results 
are presented, discussed, and compared to previous work; in 
Section 4 we consider some caveats, and in Section 5 we offer 
our conclusions. 

2. Description of the model. 

We model planetary accretion using a modified version of the 
Mercury 6 hybrid- symplectic integrator (Chambers 1999), run 
as an A'^ -I- A^' simulation, where we have A'^ protoplanets em- 
bedded in a disk of A^' "super-planetesimals" - particles with 
masses a tenth that of the initial masses of protoplanets which 
act as statistical tracers of a much larger number of real plan- 
etesimals and are capable of exerting dynamical friction on 
larger bodies (e.g. Thommes et al. 2003). The protoplanets (and 
the giant when it is introduced) interact gravitationally and can 
accrete and merge inelastically with all the other bodies in the 
simulation, whereas the super-planetesimal population is non- 
self-interacting. Super-planetesimals however are subject to a 
drag force from their motion relative to the nebular gas that is 
equivalent in its dynamical action to the gas drag that would be 
experienced by a single 10 km radius planetesimal. A detailed 
outline of these aspects of our model is given in Fogg & Nelson 
(2005). 

We model the evolution of the nebular gas using a 1-D vis- 
cous disk model that is linked to the N-body code with a time- 
step synchronization routine. This model solves numerically a 
viscous gas disk diffusion equation that is modified to include 
the tidal torques exerted by an embedded giant planet (Lin & 
Papaloizou 1986; Takeuchi et al. 1996) and its implementation 
is fuUy described in Fogg & Nelson (2007). The gas responds 



under the influence of this algorithm by depleting over time via 
viscous accretion onto the central star; by opening up an an- 
nular gap centered on the giant planet's orbit; and by forming 
a partial inner cavity due to dissipation of propagating spiral 
waves excited by the giant planet. The back reaction of these 
effects on the giant planet is resolved as torques which drive 
type n migration in a fully self-consistent manner. 

2.1. Type I migration forces. 

We have added type 1 migration forces to our model, which are 
applied to protoplanets only, by implementing the approach of 
Cresswell & Nelson (2006) who adopted the migration time 
prescription of Tanaka et al. (2002) and eccentricity damping 
time prescription of Tanaka & Ward (2004) and modified them 
with factors derived by Papaloizou & Larwood (2000) to de- 
scribe evolution in the case of large eccentricity. 
Their formula for the type I migration time is: 




where M, is the mass of the central star, 2g is the gas surface 
density in the vicinity of the planet, h is the gas scale height, 
and m, a, r, e and Q. are the planet's mass, semi-major axis, dis- 
tance from the star, orbital eccentricity and orbital frequency 
respectively. The factor /3 in the first term is the gas disk sur- 
face density profile index (2g oc r"^), which we take to be fixed 
at its initial value of fi - 1.5, even though fi falls in value dur- 
ing the evolution of the inner disk as gas drains onto the star. 
In practise, this simplification makes only a ~ 15% difference 
to t„ at later times and avoids the need for additional, time- 
consuming, measurement of yS at every protoplanetary location 
each gas disk time step. Important behavioral features to note 
from Eq. 1 are that type I migration speeds up with an increase 
in planetary mass and slows down with a decrease in gas den- 
sity or increased eccentricity. When e > 1.1 hjr, inward mi- 
gration halts as tm becomes negative and only resumes when 
eccentricity is damped to lower values. We note that under cer- 
tain circumstances, such as there being a surface density jump, 
or an optically thick disk, or MHD turbulence type I migration 
may be substantially modified or reversed (Masset et al. 2006; 
Paardekooper & MeUema 2006; Papaloizou & Nelson 2005; 
Nelson 2005). Consideration of these effects goes beyond the 
scope of this paper. 

The control of ^g{r, f) by our viscous disk algorithm 
causes fm(r, t) to behave in an interesting way. Since tm oc 
Eg'a"^(/!/r)^£l"' and Sg oc r"^, and as we use hjr - 
0.047(r/AU)'^'* AU, we find f„ oc and r oc At early 
times and in our outer disk, y8 ^ 1.5 and r oc r~^l^, so migra- 
tion is faster for interior objects of a given mass. However, the 
viscous draining of gas onto the central star rapidly results in 
much of the nebula declining to a shallower profile index of 
P - I, falling further to yS ^ 0.75 interior to 1 AU (see Fig. 1 in 
Sect. 2.2); in this situation r oc r'^'*, so relative migration rates 
are slower for interior objects. 
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Cresswell and Nelson's formula for eccentricity damping 
time is: 

where =0.1 was chosen as a normalization factor to get 
into good agreement with values measured from hydrodynamic 
simulations. We additionally assume that the timescale for in- 
clination damping is the same as that for eccentricity damping, 
i.e. ti = tg. As with the behavior for inward migration, eccen- 
tricity damping is also stronger with increasing planetary mass 
and weaker with decreasing gas density or higher eccentricity. 

Given tm and t^, type I migration forces are then applied via 
the following accelerations to each protoplanet: 

am = (3) 

a, = -2^1^ -2^, (4) 

where v is the protoplanet' s velocity vector and k is a unit vec- 
tor in the vertical direction. 



Table 1. Data describing initial solids disk set-up 





Rocky Zone 


Icy Zone 


Total 




0.4-2.7 AU 


2.7-t.O AU 


0.4-4.0 AU 




9.99 Me 


14.8 Me 


24.8 Me 


'^proto 


0.025 Me 


0.1 Me 




N 


66 


9 


75 




0.0025 Me 


0.01 Me 




N' 


3336 


1392 


4278 


./proto 


0.17 


0.06 


0.1 




100 



2.2. Initial conditions and running of ttie simulations. 

As an initial state, we assume the presence of a protoplanetary 
disk about a 1 Mg star conforming in properties and structure 
to the minimum mass solar nebula model of Hayashi (1981), 
but scaled up by a factor of three. We assign a nominal age 
of this disk to be ~ 0.5 Myr (this being the t = start time 
for the simulations) by which point the evolution of the solid 
components of the inner disk is considered to have reached the 
oligarchic growth stage described by Kokubo & Ida (2000). 
This inner solids disk extends initially between 0.4-4.0 AU and 
has a snowline at 2.7 AU beyond which its surface density is 
boosted by a factor of 4.2 to account for ice condensation. The 
generation of the N-body components of this initial solids disk 
assumes initial protoplanetary masses of 0.025 and 0.1 M® in- 
terior and exterior to the snowline respectively, spaced approx- 
imately 8 mutual Hill radii apart, with the remainder of the ma- 
terial inventory consisting of super-planetesimals with a fixed 
mass of 10% of the initial masses of the local protoplanets. 
Further details are described in Fogg & Nelson (2005) and rel- 
evant data are shown in Table 1 which gives, for zones interior 
and exterior to the snowline, values for the total mass of solid 
material Msoud, the number and mass of protoplanets A'^ and 
Wproto> and the number and mass of super-planetesimals A^' and 
Ws-pi- The parameter /proto> at the foot of Table 1, is the mass 
fraction of the solids disk contained in protoplanets and we use 
this here as a rough measure of the evolution of the disk, tak- 
ing /proto = 0.5 to denote the transition between oligarchic and 
giant impact growth regimes. We model the gas component of 
the disk from 0.025-33 AU with an initial surface density pro- 
file of Sg (X r"'-^ which gives an initial mass of 0.0398 Mq 
s 42 Mj. The disk alpha viscosity is a = 2 x 10"^ which gives 
a viscous evolution time at 5 AU - 120000 years and a mass 
depletion e-folding time ^ 580000 years (see Fogg & Nelson 
(2007) for details). 
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0.1 1 10 

r(AU) 

Fig. 1. Evolution of the gas surface density within the inner 10 AU of 
our simulated disk. The upper line is the initial r'^-^ Sg-profile for a 
3 xMMSN disk. The lower curves, in descending order, are the profiles 
at 0.1, 1.0 and 1.5 Myr respectively. 

From ; = we start by running this combined N-body and 
viscous disk model in the absence of the giant planet in order 
to mature it to different ages. This is done with a symplectic 
time step of 8 days and a simulation inner edge of 0.1 AU, in- 
terior to which any solid material is eliminated. Our previous 
simulations have shown the importance of this maturation as 
the relative strengths of dissipative forces (gas drag, dynamical 
friction, and the newly included type I migration) all vary with 
time, as the gas density falls and protoplanets grow at the ex- 
pense of planetesimals. We generate one "early" and two "late" 
scenarios by ageing the system for 0.1, 1.0 and 1 .5 Myr that are 
the equivalent of Scenarios I, IV and V from Fogg & Nelson 
(2007) and to facilitate comparison with this previous work we 
keep this nomenclature here. 

The initial gas surface density and its profile at these three 
ages are plotted in Fig. 1 and numerical data for the evolved 
disks, including the remaining mass of gas (Mg^J and solid ma- 
terial (Msoiid), the maximum protoplanetary mass (mmax)> the 
number of surviving particles (N & N'), and the protoplanet 
mass fraction (/p,-oto) are given in Table 2. Large reductions in 
gas density are illustrated, especially in the inner regions of the 
system. Particle numbers fall and mmax and /proto rise with time 
as expected and by the time of the latest scenario at f = 1 .5 Myr, 
some solid material has been lost within 0.1 AU. It might be 
thought that this loss (just ~ 5% of our solids inventory) is 
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Table 2. Disk data: after 0.1, 1.0 and 1.5 Myr of evolution 



Time (Myr) 


0.1 


1.0 


1.5 


Scenario ID 


I 


IV 


V 


Mg,, (Mj) 


30 


7 


3 


Msoiid (Me) 


24.8 


24.8 


23.5 


Wmax (Me) 


0.55 


1.07 


1.83 


N 


51 


24 


16 


N' 


4087 


2171 


1585 


.^oto 


0.20 


0.47 


0.51 



considerably less than would be expected given our inclusion 
of type I migration forces, however at early times the protoplan- 
etary masses are too small to be greatly effected by inward type 
I drift and by the time they have grown an order of magnitude to 
~ Me, inward drift remains slow as the gas density in the inner 
system has fallen by an even larger factor. The state of the solid 
components of the disk at the starting point of Scenario I shows 
little difference to its equivalent case with no type I migration in 
Fogg & Nelson (2007); however by the time the system is aged 
to the starting point of Scenarios IV & V, substantial inward 
movement of the protoplanetary population has happened, de- 
pleting the outer disk and crowding the inner regions with pro- 
toplanets (see Fig. 2). The adequacy of an 8 day time-step for 
our late maturation runs is questionable as significant material 
has moved inside 0.4 AU at f > 1 Myr. However, protoplane- 
tary growth in the inner part of our model disk occurs rapidly 
and mostly at distances of > 0.4 AU and tests we conducted 
to examine this issue showed that disk forces prevent strongly 
divergent behavior as a function of step size. The mass distri- 
bution of our matured disks interior to 0.4 AU is therefore not 
significantly affected by our choice of step size and, once the 
giant planet is inserted at the start of each scenario, orbital evo- 
lution is resolved with a shorter time-step of one tenth the or- 
bital period of the innermost object. As well as type I migration 
generating edge effects at the inner edge of our solids disk, it 
is likely that it would drive in other objects into our simulation 
region from beyond 4 AU over this time: we have not included 
this possibility here but connment on it further in Sect. 4. 

The three aged disks summarized in Table 2 are used as the 
basis for the type II giant planet migration scenarios studied 
here. In each case, a giant planet of 0.5 Mj is inserted at 5 AU 
after removing 0.4 Mj of gas from the disk between 3-7 AU 
and the inner boundary of the simulation is reset to 0.014 AU, 
which is approximately the radius of a T-Tauri star (Bertout 
1989). The giant planet then proceeds to clear an annular gap in 
the gas and migrates inward in step with its viscous evolution. 
The simulations are halted when the giant reaches hot- Jupiter 
territory at 0. 1 AU: in Scenario I, this takes ~ 1 14 000 years, in 
close agreement with the viscous disk evolution time at 5 AU, 
but in Scenarios IV and V the process takes longer (~ 152500 
and 220 000 years respectively) because at their late times of 
starting most of the gas has been lost (see Table 2) and the re- 
mainder is less effective at driving migration. The symplectic 
time-step for these runs was set to one tenth the orbital period 
of the innermost object which was achieved by dividing each 
simulation into a set of sequential sub-runs with the time- step 



2.0- 
1.8- 
1.6- 
1.4- 
1.2- 

i" 1.0- 

w 

CO 

J 0.8- 
0.6- 
0.4- 
0.2- 
0.0- 



o t = 0.1 Myr 
o t = 1 .0 Myr 
• t = 1.5 Myr 



» o 



• •• • 



•o 



o o 



. o.cn , 



1.5 2.0 2.5 

a(AU) 



Fig. 2. Protoplanetary masses vs. semi-major axis obtained after evo- 
lution to the starting points of Scenarios 1, IV and V. Substantial in- 
ward drift is visible at late times due to type I migration. 



adjusted appropriately at each re-start. At late times in each 
scenario when material had drifted interior to 0.1 AU, this en- 
tailed working with time-steps of less than one day with overall 
run times of several months as a consequence. 

3. Results. 

We describe first the behavior of one of our runs in detail and 
then proceed to discuss the effect of adding type I migration to 
our model by contrasting the results of our three scenarios pre- 
sented here with their counterparts in Fogg & Nelson (2007). 

3.1. An account of Scenario IV 

Scenario IV is one of our "late" scenarios, taking place 1 Myr 
after our model start time at a system age of ~ 1.5 Myr. By this 
point, and after removing the gas required for the giant planet's 
envelope, the amount of gas remaining in the disk has fallen to 
just 16% of its initial value with the result that the giant takes 
152000 years to migrate inward to 0.1 AU, longer than would 
be predicted from the viscous evolution time. Four snapshots 
of the evolution of Scenario IV are illustrated in Figs. 3-6 
showing the mass, incUnation and eccentricity of objects, and 
the gas surface density vs. semi-major axis. The original prove- 
nance of the protoplanets (interior or exterior to the snowline) is 
denoted by the shading of its symbol as described in the caption 
to Fig. 3. In the case of a merger between rocky and icy proto- 
planets, this shading is determined by that of the most massive 
of the pair 

An early stage in the run at 20 000 years after the introduc- 
tion of the giant planet is shown in Fig. 3. The giant has moved 
in to 4.29 AU and has caused significant excitation of the outer 
solids disk material at first order resonances. Noticeable also 
is that the outer disk beyond ~ 2 AU has been largely cleared 
of protoplanets due to previous type I migration. Only two low 
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Fig. 3. Scenario IV at 20000 years after 
the start of giant planet migration, show- 
ing the mass, inclination and eccentricity of 
objects. Small black dots represent super- 
planetesimals; white filled circles are rocky 
protoplanets; grey filled circles are icy pro- 
toplanets and the large highlighted grey 
filled circle is the giant. Objects plotted be- 
tween the dotted lines in the upper panel 
have orbits that intersect the orbit of the gi- 
ant. The location of the 2: 1, 3:2 and 4:3 reso- 
nances with the giant are indicated. Gas sur- 
face density is read on the right hand axis of 
the lower panel, the upper grey curve being 
the unevolved profile att = and the lower 
black curve being the current profile. 
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Fig. 4. Scenario IV at 80 000 years after the 
start of giant planet migration. The giant has 
now moved inward to 2.04 AU and is be- 
ginning to enter the inner zone that has be- 
come packed with protoplanets due to pre- 
vious type I migration. 



mass protoplanets remain in this region, one of them having 
recently become caught at the interior 3:2 resonance with the 
giant planet and the other having already been scattered into an 
external orbit. 

An intermediate stage in Scenario IV at 80000 years is 
shown in Fig. 4 where the giant planet has migrated inward to 
2.04 AU and is at the point of entering the zone that is crowded 
with protoplanets at a < 2 AU. Most of the planetesimals orig- 
inating in the outer disk have now been shepherded interior to 
2 AU, greatly increasing the surface density of solids caught 
between the 2:1 and 4:3 resonances. Two protoplanets have 
been captured at the 3:2 resonance and are in the process of 



having their orbits excited to an eccentricity where they will 
cross the giant planet's orbit and become vulnerable to scatter- 
ing or accretion. Type I eccentricity damping is therefore not 
always strong enough to counteract the resonant pumping ex- 
erted by the giant planet as is shown by the additional 0.46 M® 
protoplanet that has been expelled into the growing scattered 
disk. One reason for this is that /<, increases with the cube of 
eccentricity (Eq. 2) greatly reducing the resistance of type I 
eccentricity damping to strong perturbations. It is also evident 
by this stage that the giant planet is easily capable of catch- 
ing up and overtaking the protoplanets which have only moved 
inward slightly in the meantime. For example, the type I migra- 
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Fig. 5. Scenario IV at 120000 years after 
the start of giant planet migration. The gi- 
ant has now moved inward to 0.72 AU. Six 
more protoplanets and a substantial number 
of planetesimals have been scattered into the 
external disk. 





km 



Fig. 6. The end point of Scenario IV at 
152 500 years after the start of giant planet 
migration. The giant is at 0.1 AU. About 
half of the original solids disk mass survives 
in external orbits and seven protoplanets are 
found between 0.7 - 3 AU. A hot-Earth re- 
mains interior to the giant planet in an ec- 
centric and incUned orbit. 



tion time given by Eq. 1 for a 0.5 protoplanet at 2 AU, with 
e = 0, h/r = 0.056 and Eg = 150 g cm"-^ (a twelve-fold reduc- 
tion in gas surface density at this location from that at f = 0) is 
~ 6 Myr, comfortably larger than the type II migration time 
of the giant. 

Fig. 5 shows the system 120 000 years after the introduction 
of the giant planet which has now reached 0.72 AU. An exter- 
nal disk is building up. Eight more protoplanets have had close 
encounters with the giant planet, resulting in two of them being 
accreted and six undergoing external scattering. The orbits of 
the most massive and least eccentric of these bodies are seen 
to be damping fast, as would be expected from the influence 



of Eq. 2. Some of this damping might also be due to dynami- 
cal friction as a substantial number of super-planetesimals have 
also been scattered into the same vicinity. Most of the solids 
mass however remains in the shepherded disk fraction, interior 
to the giant planet. Gas drag and type I damping forces in this 
region however are now two orders of magnitude weaker than 
a nominal 3 x MMSN model as gas accretion onto the central 
star has reduced Eg here by a factor of ~ 25 since t = Q and 
dissipation of spiral waves generated by the approaching giant 
have recently reduced Zg by an additional factor of ~ 4. 

Scenario IV ends at 152500 years after the introduction of 
the giant planet when it reaches 0.1 AU and the system con- 
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Table 3. Fate of the disk mass at the end of each run: results for scenarios without (-») and with (—>) type I migration. 



Scenario 


I^* 


I^ 


IV^ 


IV^ 






Total Initial Solids (Mg,) 


24.81 


24.81 


24.81 


24.81 


24.81 


23.47 


Total Surviving Solids (Me) 


16.60 (67%) 


18.70 (75%) 


21.23 (86%) 


15.13 (61%) 


20.22 (81%) 


18.23 (78%) 


Interior Surviving Solids (M^) 


0.88 (4%) 


3.44 (14%) 


0.84 (3%) 


1.84 (7%) 


0.31 (1%) 


6.23 (27%) 


^» .^roto 


0,0 


1, 0.99 


0,0 


1, 0.96 


0,0 


1, 1 


Exterior Surviving Solids (Me) 


15.72 (63%) 


15.27 (62%) 


20.39 (82%) 


13.29 (54%) 


19.90 (80%) 


12.00 (51%) 


A',/proto 


39, 0.27 


34, 0.22 


31,0.63 


9, 0.34 


23, 0.66 


7, 0.24 


Accreted by Star (Me) 


0.01 (0.04%) 


0.0 (0%) 


0.0 (0%) 


0.0 (0%) 


0.0 (0%) 


0.01 (0.04%) 


Accreted by Giant (M®) 


8.20 (33%) 


6.01 (24%) 


3.41 (14%) 


9.68 (39%) 


4.59 (19%) 


5.23 (22%) 


Ejected (M®) 


0.00 (0%) 


0.1 (0.4%) 


0.17 (1%) 


0.0 (0%) 


0.0 (0%) 


0.0 (0%) 



28- 
27- 
26- 
25- 
24: 
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Fig. 7. Scenario rv interior to 3 AU at the point where the giant planet 
reaches 0.1 AU. Inclination is plotted vs semi-major axis. Super plan- 
etesimals are the grey dots and protoplanets the open circles. Planets 
also have marked their masses in M® and their radial excursions due 
to orbital eccentricity. 



figuration at this point is shown in Fig. 6. The final 0.6 AU 
of the giant's migration has had a destructive effect on the 
densely packed inner system material. Of the twelve proto- 
planets remaining interior to the giant at the point shown in 
Fig. 5, there are only two remaining: a 0.79 planet scattered 
to 0.8 AU, and a 1.77 hot-Earth candidate at 0.062 AU, 
close to the 2:1 resonance, with an eccentric, inclined orbit 
(e = 0.46, i = 25.3°). The other ten, which were poorly damped 
due to low gas densities, and compacted into an ever shrinking 
volume, underwent a series of violent mutual impacts followed 
by eventual accretion by the giant planet. A close encounter 
with the last of these hot-Earths to hit the giant was responsible 
for the high inclination of the orbit of the one survivor which, 
at i = 25.3°, is too low for the Kozai mechanism to play a 
role in subsequent orbital evolution. Despite these losses how- 
ever, about half the mass of the original solids disk survives in 
the scattered external disk and most of this at radial distances 
0.7 < r < 3 AU (see the magnified illustration in Fig. 7). Seven 
substantial protoplanets ranging in mass from 0.26 - 0.82 
now occupy this region along with a large quantity of yet to 



be accreted planetesimals. Type I eccentricity damping has al- 
ready circularized the orbits of four of the protoplanets which 
might be expected to undergo some degree of renewed inward 
type I migration before complete gas dispersal followed by a fi- 
nal phase of planetesimal clear- up and giant impacts. Allowing 
for this, the configuration shown in Fig. 7 suggests a good 
chance that the final planetary arrangement here could include 
a planet in the system's habitable zone and we return to this 
issue in Sect. 3.3. 



3.2. Comparison with simulations not including type I 
migration. 

Data showing the fate of the solids disk mass at the end of the 
three scenarios presented here that include type 1 migration are 
given in Table 3 in the columns headed with a — » symbol and 
are compared against their counterparts from Fogg & Nelson 
(2007) in the columns headed with a -f* symbol. Each column 
gives the total mass of solid material in the disk at the start 
of giant planet migration; the mass and the percentage of the 
initial disk surviving at the end of the migration; and the mass 
surviving interior and exterior to the giant planet's final orbit 
at 0.1 AU, including the number of remaining protoplanets and 
the protoplanetary mass fraction in each respective partition. 
The bottom three rows give the data for the mass loss channels 
which are accretion by the central star, accretion by the giant 
planet and ejection into a hyperbolic orbit. 

These data show that a large fraction of the disk mass sur- 
vives the passage of the giant planet in each case: inclusion 
of type I migration at the nominal rate does not result in the in- 
ner system being cleared of planet-forming material. In the two 
versions of Scenario I slightly more mass remains with type I 
migration operating than without. However, this is largely due 
to the presence of a 3.42 M® hot-Earth inside 0.1 AU captured 
in the 3:2 resonance with the giant planet with an eccentric- 
ity of 0.29. Generally the data for Scenarios 1^ and 1^ are 
similar, suggesting that type I migration only has a small influ- 
ence on the outcome of giant planet migration through a young 
inner system disk where protoplanetary masses remain small. 
In our more mature Scenarios IV & V we do find less solids 
mass surviving when type I migration is included, but this loss 
is not catastrophic: 61% of the solids disk remains at the end 
of Scenario IV (see Figs. 6 & 7) and 78% remains at the end 
of Scenario V. Almost exclusively, the mass that is lost is ac- 



MJ. Fogg & R.P. Nelson: Effect of type I migration on formation of terrestrial planets in hot- Jupiter systems 



9 



creted by the giant planet. The increase in efficiency of this loss 
channel results from the fact that, due to previous inward type I 
migration, protoplanetary encounters with the giant planet tend 
to occur at smaller orbital radii where the probability of accre- 
tion versus scattering is relatively enhanced. In Scenario V this 
additional mass loss is minor as much of the inner system mate- 
rial had by this late time accreted into substantial protoplanets 
which were able to migrate ahead of the giant inside 0. 1 AU: a 
situation that can be seen developing in Fig. 2. These eventu- 
ally co-accreted to form a 6.23 hot-Earth that was captured 
into the 2:1 resonance with the giant planet near the end of the 
simulation exciting its orbit to an eccentricity of 0.41. Future 
accretion of this body by the giant planet is a possibility and 
we examine the issue of the long-term survival of our hot-Earth 
planets in Sect. 3.4. 

Across the board in Table 3, the consistent trends that 
emerge when type I migration is added to our model are: 

1. Shepherding is modestly enhanced. This can be seen from 
the larger percentage of the original solids disk mass that 
is found in the interior surviving sohds fraction, which in- 
cludes in each case a single hot-Earth type planet with an 
eccentric orbit captured at a first order resonance with the 
giant planet. Such objects were seen to form, but not to 
survive, in our analogous runs in Fogg & Nelson (2007): 
their typical fate being accretion by the giant planet. The 
hot-Earths generated here have narrowly escaped this fate 
themselves and their long-term survival would depend on 
the cessation of any further inward migration of the giant 
and the effectiveness of any continuing type I eccentricity 
damping exerted by the tenuous gas (reduced by over two 
orders of magnitude from a 2^ oc r"'-^ model) remaining 
in the planet's vicinity. We do not recover the well-damped 
hot-Neptune or hot-Earth type planets generated in Fogg & 
Nelson (2005) where we assumed a steady state gas disk 
with a fixed oc r"^^ surface density profile. In this case, 
the persistently high gas densities caused more rapid or- 
bital decay of planetesimals, strong dynamical friction and 
high rates of collisional damping which were sufficient to 
drag several protoplanets into low orbits where a rapid final 
phase of accretion was often sufficient to free the remain- 
ing body from resonant locking with the giant planet. Here, 
our evolving and decaying gas disk weakens both dynam- 
ical friction and type I eccentricity damping; whilst type 
I orbital decay facilitates inward movement of protoplan- 
ets, at late times it also brings about a separation of pro- 
toplanets from the bulk of planetesimals, further reducing 
the influence of dynamical friction on the innermost bod- 
ies. The overall dissipation in our present model therefore 
only results in a maximum of 27% of the original solids 
disk partitioned into the inner remnant. 

2. Scattering is modestly reduced and strongly selective. Our 
simulations that include type I migration produce less mas- 
sive scattered disks, but not dramatically so. A stronger ef- 
fect is that the likelihood of external scattering is now bi- 
assed in favor of planetesimals, generating immature exte- 
rior solids disks with fewer protoplanets. This can be seen 
in Fig. 8 where the black bars show the mass of the scat- 




Fig. 8. Mass fractions of scattered disks generated in the present work 

compared to those of analogous scenarios not including type I migra- 
tion and before and after protoplanet mass fractions. 



tered disk resulting from runs with no type I migration and 
the grey bars show those resulting from the present work. 
The superimposed lozenge symbols show the disk proto- 
planet mass fraction (fproto) at the start of each scenario and 
the circle symbols give the value of fproio at the end. In late 
scenarios the addition of type I migration to our model has 
resulted in the reduction of the mass of the outer disk rem- 
nant by ~ 30% which still leaves ~ 50% of the mass of the 
original ending up in the scattered disk. These disks how- 
ever are less mature in the sense that they contain a greater 
fraction of small bodies. In Fogg & Nelson (2007) it was 
noted that one effect of giant planet migration through a 
solids disk was to increase fproto in the external renmant by 
selectively scattering protoplanets. Fig. 8 shows that in late 
scenarios, this trend is reversed here as type I migration, 
especially at late times, causes protoplanets to be prefer- 
entially shepherded: encounters with the giant tend to oc- 
cur closer in so that fewer of them escape into the exter- 
nal renmant allowing the scattered planetesimal population 
to dominate. The external disks generated here however, 
though of reduced mass, should still be capable of support- 
ing future planet formation. The protoplanets they contain 
can be of substantial mass and are typically emplaced at 
smaller semi-major axes and the more abundant planetes- 
imal population and residual type I migration forces pro- 
vide for stronger damping of protoplanets after scattering 
has occurred. 



3.3. Post-migration terrestrial planet formation 

To model the completion of planet formation in our scattered 

disks would require simulating a further ~ 100 Myr of accre- 
tion and is beyond the scope of this paper. Shorter range in- 
tegrations however are still of considerable interest as type I 
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Fig. 9. Details of the inner 3 AU at 
the end of Scenario IV (top panel) and 
1.3475 Myr later at the point of gas dis- 
persal. Eccentricity is plotted against semi- 
major axis with orbital excursions due 
to eccentricity illustrated with error bars. 
Planetary masses are labeled in units of Mg,. 
Eg is read from the right hand axes: the up- 
per curve being the profile at f = (sys- 
tem age = 0.5 Myr) and the lower curve 
being the present profile. The extent of the 
maximum greenhouse habitable zone (0.8 - 
1.7 AU) is shown by the horizontal bar. 



migration forces only decline to zero when aU the remaining 
local gas has gone. In this section, we continue the evolution of 
Scenario IV to a system age of 3 Myr, requiring an additional 
1 .3475 Myr of simulated run time, at which point we assume all 
residual gas is rapidly lost to photoevaporation. Over this pe- 
riod the mass of the gas disk falls from ~ 5.5 Mj to ~ 0.55 Mj. 
In order to evade the severe limitations of our viscous disk al- 
gorithm time-step, which adapts to a tiny value with the giant 
planet at ~ 0.1 AU, we turn off the type 11 migration forces 
effecting the hot- Jupiter and replace our viscous disk with an 
exponentially declining disk that maintains the previous sur- 
face density profile and declines with a mass e-folding time of 
581 514 years. We evolve the system henceforth with a 2 day 
N-body time-step. This set-up assumes an unspecified migra- 
tion halting mechanism on the part of the hot- Jupiter that is not 
fully consistent with our model and we comment on the conse- 
quences of this in Sect. 4. 
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Details of the evolution of the inner 3 AU of the sys- 
tem during this terminal phase of gas depletion are shown in 
Fig. 9 where the top panel shows the state of play at the end 
of Scenario IV and the bottom panel shows it at 1 .3475 Myr 
later at a system age of 3 Myr. During this period we find that 
dynamical friction and residual type I migration forces rapidly 
damp and circularize the orbits of protoplanets with high eccen- 
tricities. Gas drag is also sufficient to maintain the eccentricities 
of planetesimals at a low value. The four well damped proto- 
planets between 1 .25 - 2.5 AU in the upper panel of Fig. 9 grow 
by a few percent from the surrounding planetesimal field and 
migrate inward an additional ~ 0.6 AU. By the time we remove 
the remaining gas in the bottom panel, two more planets have 
moved into the maximum greenhouse habitable zone (Kasting 
et al. 1993) giving a total of three planets sited between 0.8 - 
1.7 AU. The migration of these four protoplanets shows some 
of the convoying behavior described in McNeil et al. (2005): 



Fig. 10. Resonant angles for the 5:2 mean motion resonance between 
the giant planet and the innermost external protoplanet. Resonant cap- 
ture is indicated ~ 125 000 years before the stage shown in the bottom 
panel of Fig. 9. 



convergent migration can cause protoplanets to bunch together 
and lock into a stable arrangement of mean motion resonances 
whereupon the pattern then migrates inward at some average 
drift rate. In the bottom panel of Fig. 9 the 0.76 M® planet and 
its 0.89 Me neighbour are caught at their 7:6 resonance and the 
outermost of this pair is close to the 4:3 resonance with its outer 
0.74 M® neighbour. The innermost external protoplanet, hav- 
ing only recently been scattered, starts with mproto - 0.78 Me, 
a ^ 0.8 AUande ^ 0.6 but damps rapidly and migrates inward, 
growing to 0.94 M® and eventually halting at a = 0.19 AU at 



MJ. Fogg & R.P. Nelson: Effect of type I migration on formation of terrestrial planets in hot- Jupiter systems 



11 



0.7- 

0.6- 
0.5- 
0.4- 
<0 0.3- 
0.2- 
0.1 ■ 
0.0- 



System age = 4.0 Myr 



"3 'i. ■ <yc 



0.7 
0.6 
0.5- 
0.4- 
0.3- 
0.2 



0.0 



0.5 



1.0 



1.5 



2.0 



2.5 



3.0 



System age = 6.0 Myr 



Maximum greenhouse habitable zone 



OA-® S 
0.0--^^ 



0.0 



— I — 
0.5 



4D- 



e 



1.0 



1 

1.5 

r(AU) 



2.0 



a- 1 



2.5 



3.0 



Fig. 12. Details of the inner 3 AU of the 
Scenario IV system after 3 Myr of gas-free 
d3mamics. The top and bottom panels show 
configurations at system ages of 4 and 6 M5T 
respectively. 
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Fig. 11. Material composition of protoplanets in the extended Scenario 
IV at the stage shown in the bottom panel of Fig. 9. 



the 2:5 resonance with the hot- Jupiter (see Fig. 10). This behav- 
ior is an example of the tendency, when convergent migration 
is in play, for giant planets to collect smaller bodies at external 
mean motion resonances (Thommes 2005). 

Our model also tracks the volatile composition of proto- 
planets, recognizing three crude material provenances: mate- 
rial originating at a < 2 AU is assumed to be dry rock; from 
between 2 AU < a < 2.7 AU its water content is assumed 
to be similar to that of chondritic asteroids; and from beyond 
a > 2.7 AU it is assumed to be water-rich trans-snowline mate- 
rial. Previous models have demonstrated how a migrating gi- 
ant planet drives volatile-rich material into its iimer system. 



boosting the potential water content of any planets that form 
there (Raymond et al. 2006; Fogg & Nelson 2007; Mandell 
et al. 2007). This trend is accentuated when type I migration 
forces are included as volatile-rich protoplanets can now mi- 
grate directly inward. In Fig. 11 we show the composition of 
the protoplanets in the extended Scenario IV at a system age 
of 3 Myr. All the bodies have a substantial endowment of outer 
solids disk material especially those at a > 0.8 AU. The largest 
three of these planets, the resonant trio situated between 1.3 - 
1.9 AU, all originated from beyond the snowline at a > 2.7 AU 
and accreted copious volatile-rich material during their inward 
journey. Ocean worlds are a probable outcome of this kind of 
evolution (Kuchner 2003; Leger et al. 2004). 

We have evolved this system, in the absence of gas, for a 
further 3 Myr to a system age of 6 Myr and two configura- 
tions over this period are shown in Fig. 12. The cessation of 
type I migration forces causes a shift to a more oligarchic style 
of growth (Kokubo & Ida 1998). Orbital repulsion breaks up 
the resonant convoy showing in Fig. 9 and the entire proto- 
planet chain dynamically expands. The cessation of gas drag 
results in planetesimal orbits becoming much more excited in 
the process of exerting dynamical friction. Although this re- 
duces protoplanetary accretion cross-sections, protoplanets no 
longer shepherd gaps in the solids disk and continuing growth 
is in progress. 

The final state of the system shown in the bottom panel of 
Fig. 12 cannot be predicted in detail. We speculate that plan- 
etesimals will eventually be accreted or cleared, with some 
small fraction of them possibly surviving in the prominent gap 
in the planetary distribution at ~ 0.5 AU. Future mergers of 
protoplanets via giant impacts are less certain as they all lie in 
near-circular and well-spaced orbits: the closest being the two 
bodies between 1-1.5 AU which are nonetheless separated by 
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Fig. 13. Further evolution of the orbits of 
the interior planets resulting from the three 
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~ 19 mutual Hill radii. We expect therefore that the final con- 
figuration of the Scenario IV system will contain either one or 
two Earth-mass planets in its habitable zone and longer range 
integrations of this system are underway. 

3.4. Hot-Earth survival 

Having used a 2 day time-step, the additional simulation repre- 
sented by Fig. 12 does not adequately resolve the orbital evo- 
lution of the hot-Earth planet at 0.062 AU. Single hot-Earths, 
captured at interior first order mean motion resonances resulted 
from all our scenarios presented here, but their eccentric orbits 
with apastra close to the radial distance of the giant planet raise 
questions over long term survival. We have examined this issue 
by running the orbits of our hot-Earth candidates for a further 
10000 years (~ 10^ orbits at 0.05 AU), with a half day time- 
step, in the absence of gas and including the giant planet and 
any remaining debris interior to 0. 1 AU. This assumes the exis- 
tence of an interior gaseous cavity which is not fuUy consistent 
with our scenario end points and we comment on this further 
in Sect. 4. 

The orbital evolution resulting from these extended runs is 
illustrated in Fig. 13 and their final configuration is shown in 
Fig. 14. All the interior planets remain. The least convincing 
as a longer-term survivor is the 3.43 Me hot-Earth present in 
Scenario I: it accreted very little additional debris (which was 
accreted by the giant instead) and continues to orbit at the 3:2 
resonance with an eccentricity of ^ 0.3, bringing it very close to 
crossing the giant's orbit at apastron. Close encounters have so 
far been avoided by the orbital phasing resulting from the res- 
onance, but we have not investigated how long this state could 
persist. In the case of Scenario IV, we start with a 1 .77 hot- 
Earth at the 2: 1 resonance in an excited orbit with e = 0.46 and 
i = 25°, but this planet succeeds in accreting an extra 0.03 
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Fig. 14. Hot-Earth - hot- Jupiter systems resulting from an additional 
10 000 years of simulation. Error bars show the radial distance be- 
tween periastron and apastron. Numeric labels give masses in M^. 



of debris, cleaning up all of it but the one remaining super- 
planetesimal showing in Fig. 14, with the result that its eccen- 
tricity is damped to e 0.2 by the end of the run period. Its 
inclination is not damped and one accretion event increases it 
slightly to i =^ 30° but the hot-Earth and the hot- Jupiter are 
now well separated and long-term survival looks promising. In 
Scenario V there was no debris for the hot Earth to clean up 
and Fig. 13 shows that its orbit remains unchanged over 10 000 
years. This 6.23 planet remains close to the 2:1 resonance 
in a very low inclination orbit with e ^ 0.41 which at no point 
comes dangerously close to the orbit of the giant planet. This 
hot-Earth is also a possible long-term survivor, especially as 
stellar tidal forces would be expected to gradually circularize 
its orbit over much longer timescales. 
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4. Discussion. 

Any simulation of planet formation is restricted by its assump- 
tions and incomplete input physics. Below, we discuss three 
areas where future improvements are possible as well as the 
caveats they might hold for our present results. 
(0 Migration halting mechanisms. In all our simulations we 
have halted the migration of our giant planet at 0. 1 AU where 
we assume that whatever mechanism is responsible for strand- 
ing hot-Jupiters in their final orbits comes into play. This mech- 
anism is unknown and proposed explanations include the pres- 
ence of a central magnetospheric cavity in the gas which could 
decouple an intruding giant planet from the disk (Lin et al. 
1996); another is that of fortuitous gas disk dispersal where 
the eventual loss of the gas disk strands those migrating gi- 
ant planets that still remain at arbitrary orbital radii (Trilling 
et al. 1998). The former could apply to any of our scenarios, 
early or late, whilst the latter is more appropriate to late sce- 
narios when most of the gas has already dispersed. Neither of 
these mechanisms are fully consistent with our present model. 
If we were to run it further, we would expect the giant planet 
to come to rest somewhere between 0.02-0.1 AU but we have 
yet to model any of the physics of a magnetospheric cavity and 
time-step limitations that stem from our viscous disk algorithm 
make it impractical to simulate further type II migration inte- 
rior to 0. 1 AU. The only difference that stopping the hot- Jupiter 
closer to the star would make to our results is that more com- 
pact systems would be less likely to host interior hot-Earth or 
hot-Neptune type planets. Further resonant shepherding of such 
objects inward of 0.1 AU runs the risk of driving them into the 
central star (Fogg & Nelson 2005). Our results for scattered 
disk characteristics and their prospects of renewed terrestrial 
planet formation would be unaff"ected by the final fate of the 
hot- Jupiter. If fortuitous disk dispersal is assumed to be the 
halting mechanism, then this could be simulated in our model 
by including a photoevaporation algorithm which removes gas 
at some small and steady rate that accords with observations. 
This would result in a rapid dispersal of the of the remaining 
gas at late times (e.g. Clarke et al. 2001; Alexander et al. 2006) 
and a termination of type II migration. Work is ongoing to in- 
clude this mechanism in our model and to conduct simulations 
where giant planets strand at a variety of distances in a fully 
self-consistent way. 

(;;) Supply of outer disk material. In our previous models the 
radial drift of protoplanets is due to dynamical stirring and a 
degree of coupling to the slow inward movement of planetesi- 
mals via dynamical friction. In these cases it was reasonable to 
assume that our giant planet formed at 5 AU and that any sup- 
ply of extra material over our run times into the inner disk from 
beyond that distance would be insignificant. However the inclu- 
sion of type I migration makes the heaviest and most evolved 
objects more inwardly mobile. As shown in Fig. 2, our late 
matured disks are largely cleared of protoplanets outward of 
~2 AU, making one wonder if this region could have been re- 
populated in the interim by objects which migrated from be- 
yond 4 AU. This seems plausible as the giant planet's core 
must also have migrated so that we can no longer propose con- 
sistently that it accreted its core at 5 AU. We have not yet at- 



tempted to model the extent of any such addition of external 
mass to our solids disk but it is clear from our results that any 
additional protoplanets encountered by the giant between 2- 
4 AU are more likely to be scattered than accreted. Thus, any 
late supply of large bodies is likely to contribute to the pop- 
ulation of the scattered disk boosting both its mass, maturity 
(fproto), and volatile content. 

(///) Planetesimal size evolution. We have assumed a uniform 
population size of 10 km planetesimals which results in an 
identical efl'ect of gas drag on any planetesimal at a particu- 
lar orbital distance. In reality, but more challenging to adapt to 
+ A^'-type simulation, we would expect an evolving plan- 
etesimal size spectrum and a much more varied response of 
small objects to gas drag. In areas of the solids disk where plan- 
etesimals are strongly stirred, such as between the 2:1 and 4:3 
resonances with the giant planet, high energy mutual collisions 
might reduce planetesimals down to smaller chunks of rubble. 
This rubble would be subject to stronger dissipation from the 
gas resulting in it being shepherded inward of the giant more 
efficiently. In Fogg & Nelson (2007) we suggested that much of 
this material should still end up in the scattered disk as it would 
be expected to accrete very efficiently onto inner protoplanets 
(see Chambers (2006b)) that would be scattered later. This ar- 
gument is weaker when type I migration is included. Supplying 
inner planets with an enhanced feedstock would not only boost 
their growth rates, but also their inward drift. We have shown in 
Sect. 3.2 how the influence of type 1 migration results in fewer 
protoplanets surviving in external orbits so boosting inward mi- 
gration rates by enhancing mass growth would accentuate this 
trend. However, as accretion rates depend as rhproto °^ we 
would expect that the bulk of this additional growth would take 
place at late times when the giant planet has compacted the 
shepherded fraction within a small orbital radius. It seems un- 
Ukely therefore that the additional dissipation resulting from 
planetesimal fragmentation would prevent entirely the scatter- 
ing of some large bodies into the external disk. 
{iv) Initial nebular surface density profile. In common with 
many other studies of terrestrial planet formation, our models 
adopt an initial condition of a MMSN-type protoplanetary disk 
with an r~^-^ surface density profile (Weidenschilling 1977; 
Hayashi 1981). Other models however predict that a shallower 
profile (E ~ r ') may be more realistic (e.g. Davis 2005; 
Andrews & Williams 2007). We have not yet explored the ef- 
fect of running our present model with a shallower initial pro- 
file but expect that it would not drastically change our results. 
Previous models that have explored this issue (e.g. Chambers 
& Cassen 2002; Raymond et al. 2005b; McNeil et al. 2005) 
find that the architectures of their resulting terrestrial plane- 
tary systems depend only weakly on the slope of the initial sur- 
face density set up with shallower gradients tending to result in 
more substantial planets forming further from the central star. 
If a giant planet were to migrate through such a disk, we expect 
that more mass would be encountered earlier in the migration, 
resulting in a more massive scattered disk and hence even more 
favorable prospects for external terrestrial planet formation. 
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5. Conclusions. 

Previous models of a giant planet migrating through an in- 
ner system protoplanetary disk have all predicted that the pro- 
cess generates a scattered external disk massive enough to sup- 
port renewed planet formation (Fogg & Nelson 2005, 2006; 
Raymond et al. 2006; Fogg & Nelson 2007; Mandell et al. 
2007). None of these works included the possible action of 
type I migration forces exerted by the gas which would cause 
heavier bodies to drift progressively closer to the central star. 
Since this process would clearly act to counter scattering we 
have questioned if its operation could invalidate previous con- 
clusions concerning subsequent terrestrial planet formation. 
Having updated our model to include type I migration and hav- 
ing re-run three scenarios from Fogg & Nelson (2007), we con- 
clude that the extra dissipation generated is still not sufficient 
to prevent the buildup of a scattered disk ~ 50% the mass of 
the original. Although these scattered disks are of lower mass, 
they are generally less dispersed and better damped, supporting 
a more rapid renewal of planetary growth. As shown in our ex- 
tended run, whilst the gas remains, residual type I forces can act 
on scattered protoplanets to draw them back into the ~ 1 AU 
region and the habitable zone. Our model also generates plan- 
ets of several Earth masses in potentially stable orbits interior 
to the final orbit of the giant planet and is thus supportive of the 
proposals of Fogg & Nelson (2005) and Zhou et al. (2005) that 
hot-Earths and hot-Neptunes found interior to hot- Jupiters may 
have originated from material shepherded inward by a migrat- 
ing giant planet. 

Hot- Jupiter systems therefore, far from being devoid of 
other planetary bodies within their warmer regions, may be 
worthwhile targets for future investigation. 
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